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• 1—1 . 

X : 1. INTRODUCTION 

The idea of a discrete space-time has been introduced by physicists in the past in 
several different ways. Heisenberg advocated a fundamental length, Snyder has proposed 
a position operator with discrete spectrum, Ahmavaara has carried out a physical model 
on a cubic lattice, which is embedded in a finite linear space over a Galois field. 

Recently the use of discrete space-time variables in Lattice Gauge Theories have 
become very popular to construct models based on experimental date. An other succesful 
line of research deals with integrable models on discrete space-time. 

Many other hypothesis can be mentioned but they are reduced to two types: i) those 
using the space-time lattice as an artificial model, and ii) those proposing the discreteness 
of space-time as some fundamental structure of reality (Lorente 1986c). 

The first assumption is very old and of technical character and do not present 
conceptual problems. The second assumption is more recent, although is based on some 
old pressupositions, namely, the relational character of the structure of space and time in 
contradistinction to the absolute idea of space and time (Earman 1989). Both assumptions 
were defended in XVII century by Leibniz and Newton, respectively, and although the last 
one become more succesful the former remained as a consistent model (Jammer 1969). 
Acording to Leibniz "Time is the order of monads not existing simultaneously. Space is 
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the order or monads that coexist or exist simultaneously. A point changes its position 
when it changes its relations from some points to different ones" ^. This assumption leads 
to some discrete structure, although Leibniz did not carried out all its consequences. Also 
Riemann is his 1854 Inaugural Dissertation discussed the posibility of a discrete manifold 
based on some intrinsic metric and Weyl took over Riemann's ideas (Griinbaum 1977). 
The relational character of space and time have been advocated recently by R. Penrose 
(1971), D. Finkelstein and C.F. von Weiszacker (Castell 1986). We have also proposed 
some physical model based on the relational hypothesis of the space and time (Lorente 
1974, 1976, 1986a,b, 1987). 

In this paper we explore some mathematical consequences of the assumption of a 
discrete space-time in his most naive structure, namely, a hypercubic lattice. Obviously 
the mathematical tool must be the functions of discrete variables and difference equations. 
Therefore, the fundamental laws of physics have to be translated into the language of 
discrete mathematics. The objection which is usually raised against such discrete schemes 
is that they are not invariant under the Lorentz group. We try to overcome this difficulty, 
as Schild did (Schild 1949), finding all integral transformations that leave the lattice of 
any dimension invariant, and applying these transformations to the fundamental field 
equations. 

In section 2 and 3, we introduce the Cayley parametrization of the classical groups 
with some easy examples that can be enlarged to all semisimple Lie groups. In section 4 
we review the isomorphims between real forms and its explicit calculation. In section 5 
we describe the method to calculate all integral transformations that leave a hypercubic 
lattice invariant, using the results of sections 2 to 4. In section 6 and 7 we construct a 
lagrangian formalism for the Klein Gordon and Dirac field on the lattice. In section 8 
we introduce the generators of the space and time displacement and Lorentz boost on the 
lattice to check the Lorentz invariance of the model and we "integrate" these generators via 
some "Taylor expansion" . From the physical point of view these schemes have a twofold 
interpretation. Either we take the continuous limit in order to recover the continuous 
character of physical laws or we keep the formulas as they appear as a consequence of the 
hypothesis of a discrete space-time. 



2. CAYLEY'S RATIONAL PARAMETRIZATION OF SEMISIMPLE 
GROUPS 

Let <S be a semisimple Lie group of complex matrices S, which leaves invariant 
some non-degenerate bilinear form. We call a matrix S of the semisimple group S non- 
exceptional if det {E + S) 7^ 0, where E is the unit matrix. Cayley (1846) has proved that 
every non-exceptional matrix S can be expressed as follows 

S={E + X)-\E - X) = {E - X){E + X)-^ (2.1) 

where X is also a non-exceptional matrix. 

^ According to Jammer, Leibniz's Monadology was inspired by the philosopher Mai- 
monides who in his Guide for the Perplexed describe some discrete structure of space and 
time. 
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If G is the coefficient matrix of the non-degenerate bihnear form, which is left invariant 
under the group S, the non-exceptional matrices S satisfy the relation 

S^GS = G (2.2) 

and because of (2.1) the corresponding matrices X will also satisfy 

X^G + GX = (2.3) 

In order to obtain the independent parameters of the semisimple group S it is more 
convenient to work with expression (2.3), which is linear, rather than with expression 
(2.2), which is quadratic. If we diagonalize or reduce to the canonical form the coefficient 
matrix G we have a further simplification of (2.3). Taking the independent elements of the 
matrix X given by (2.3) to be the independent parameters, we obtain Cayley's rational 
parametrization of the semisimple group S. (Note that when the independent elements of 
X are complex their real and imaginary part should be taken as independent parameters.) 

In Table 1 we give the explicit conditions on the non-exceptional matrices S and X for 
all classical groups, as derived from expression (2.2) and (2.3), respectively. The notation 
means the transpose matrix and 5^ the adjoint. Also 

\E + S\= det{E + S) 

and 

where E^, Ep, Eg are the unit matrix of order n,p,q respectively. The condition on the 
matrix X gives automatically the unimodularity condition 

det{E + S) = det(E - S) (2.4) 

except in the groups SU (n+1) and SU{p, q) and therefore (2.4) imposes an extra condition 
on the parameters corresponding to these groups. 



TABLE 1. Cayley's decomposition of semisimple groups 





Conditions 


Conditions 






Group 


on S 


on X 


Unimodularity 


Parameters 


S0{2n) 


S'^S = E 


X^ + X = 




n{2n - 1) 


S0{2n+1) 


S'^S = E 


X^ + X = 




n(2n+ 1) 


SU{n + 1) 


S^S = E 


xt + X = 


\E + S\ = \E-S\ 


n(n + 2) 


Sp{2n) 


S^JS = J 


X'^J + JX^O 




n(2n + 1) 


SO{p,q) 


S^IS = I 


X^I + /X = 




l/2{p + q){p + q 


SU{p,q) 


S^IS = / 


X^I + /X = 


\E + S\ = \E-S\ 


{p + g)2 - 1 
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3. SOME EXAMPLES 



3.1. The Rotation Group, SO{3) 

From Table 1 the matrix X is antisymmetric and it can be expressed in the foUowing 



way 

^ " n 




q (3.1.1) 
-q J 

where n,p,q are independent parameters and m has been introduced for convenience. 
Using (3.1.1.) one obtains the Cayley parametrization ot the non-exceptional matrix of 
the rotation group 

s = 



+ + + q^ 

m? — 71^ — p^ + —2mn + 2pq 2mp + 2nq 

X I 2mn + 2pq m? — ii? + p^ — q^ —2mq + 2np | (3.1.2) 

—2mp + 2nq 2mq + 2np -\- — p^ — q^ 

If we define a = m + m, (3 = p — iq and then impose m? + n'^ + p'^ + q'^ = 1, the 
parametrization of the matrix S given (3.1.2) is identical with the parametrization used 
by Wigner (1959) for the 3-dimensional rotation group. The parameters a and /? used by 
him are related to the parametrization of SU{2), the covering group of S0{3), in this way 

" ^ ' \a\' + \P\' = l (3.1.3) 



p* a 

In terms of the components of the axis of rotation (01,02,03) and of the angle of 
rotation 0, the Cayley parameters have the following geometrical interpretation 

q p q - in? - p^ - q^ a ^ a\ 

— = — = — , cos0= — ^ ^ ^ ^ (3.1.4) 

oi 02 03 + + + 

Similar parametrization and corresponding properties can be obtained for the N- 
dimensional rotation groups. 



3.2. The Unitary Group SU{2) 

From Table 1 the matrix X is antihermitian and it can be expressed as 
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^=7(-;' (^-^-i) 

where a and h are real parameters, p = r + is and I has been added for convenience. 
From (2.5) and (3.2.1) one obtains 

^ _ ]_ f P + ab-\p\^ + i2lb -2lp \ . 

A\ 2lp* i^ + ab-\p\'+i2la) ^"^ 



with A = l^ -ab+ |p|^ + iP{a + b). 
The antihermiticity of X gives 

det (E + X)* = det (E - X) 

but it does not imply the unimodularity condition. (In the rotation group the antisymme- 
try of X does imply the unimodularity of S.). If we impose det S = 1, from (2.1) follows 
that 

det {E + X) = det {E - X) 
Both conditions, unitarity and unimodularity of S, give det(£^ + X)=real, or 

a + b^Ti X = (3.2.3) 

Substituting (3.2.3) in (3.2.2) we obtain the general expression for the unitary 
unimodular matrices in two dimensions 

}_ ( l'^ - a? - r'^ - - i2la -2lr - i2ls 

A V 2Zr - i2ls l^-a^-r^-s^ + i2la 

with A = P + 0^ + + s^. Obviously the matrix (3.2.4) is equivalent to (3.1.3), but uses 
different parametrization. 



(3.2.4) 



3.3. The Proper Lorentz Group SO(3.1) 

From Table 1 one obtains the traceless matrix 



X 



m 






n 


-P 


r 


— n 





Q 


s 


P 







t 


r 


s 


t 






(3.3.1) 



where n,p, q, r, s, t are real independent parameters and m has been introduced as before. 
The unimodularity of S does not impose further conditions on these parameters. From 
(2.1) one gets 



1 

A 



-n'^ -p'^ + + - -t^ + 
2mn + 2pq + 2rs — 2 At 
-2mp + 2nq + 2rt + 2As 
—2mr — 2ns + 2pt — 2Xq 

-2mn + 2pq + 2rs + 2At 
m? - 11? -\- p^ - q^ - r"^ -\- - + A^ 

2mq + 2np + 2st — 2Xr m' 
-2ms + 2nr - 2qt - 2Xp 



2mn + 2nq + 2rt — 2Xs 
-2mq + 2np + 2st + 2Ar 
+ v? -p^ -r"^ - s"^ + + X^ 
-2mt - 2pr + 2qs - 2 An 



-2mr + 2ns - 2pt - 2Xq 
—2ms — 2nr + 2qt — 2Xp 
—2mt + 2pr — 2qs — 2Xn 



w? + -\- p^ + q^ 



+ + 



+ A^J 



(3.3.2) 
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where 

m\ = nt + ps + qr 

and 

A = + + / + 5^ - - - - 

If A > 0, since det S = 1, one obtains the general expression for the non-exceptional 
matrices of the proper Lorentz group (5*44 > 0). 

If r = s = t = 0, one recovers expression (3.1.2) for the proper orthogonal group in 
3- dimensions. 

lfn = p = q = one is left with the non-exceptional matrices of the pure Lorentz 
transformations. In this case, comparison of (3.3.2) with a pure Lorentz transformation 
with velocity v along v gives (M0ller, 1952) 

Vfr. Vy Vz 2mc f ^ \ — r"^ — — 



t + r'^ + s'^ +t^' \^ c^J 777,2 + r2 + s2 + ^2 (3-3.3) 



where c is the velocity of light in vacuum. 
If we define 



a = m — t + i{n — \)^ (3 = —p — r + i {q — s) 
^ = p — r + i{q + s), 6 — m + t — i{n + X) 



(3.3.4) 



and introduce these variables in the general expression of the proper Lorentz group in 
terms of the parameters of the SL{2,C) group (Naimark, 1964a) 

r I aS — 3^ = 1 
1 S J 

we obtain the expression (3.3.2) plus the condition 

mX =nt + ps + qr \ (qq^\ 

A =m^+n^+p^ + q^-r^-s^-t^-X^^lj ^^"^-^^ 



4. ISOMORPHISM BETWEEN REAL FORMS 

According to Cartan theory, there are some real forms of simple Lie groups of low 
dimensionality which are locally isomorphic (Helgason, 1978). We describe them by the 
bijection of i?" onto a set of matrices S. 

i) SL{2, R) ^ 5*0(2, 1). Define a set of 2x2 real matrices A, by the conditions = A, 
where A'^ means transposed. The bijection of an element (xq, xi,X2) of onto a matrix 
A is the following: 

^^(xo + x, X, \ 

\ Xi Xq - X2 J 

The transformations A' = SAS'^ with S e SL{2, R), map A into itself. Since 

det A = xl-xl-xl = det A' (4.2) 
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this transformation induces the desired isomorphism. 

(ii) SL{2, C) SO{3, 1). Define A, a 2x2 complex matrix, by the condition = A, 
where means the Hermitian conjugate matrix. The bijection of {xo,xi,X2,X3) in i?^ 
onto A is given by 

' Xo+ Xs Xi + 1X2 
Xi — 1X2 — X^ 



(4.3) 



The transformation A' = SAS'^ with S e SL{2, C) maps A into itself, as it is well 
known (Gel'fand et al., 1963). Since 



ddj-Zl — 37q X-y — d-Gt-/! 

this transformation induces the mentioned isomorphism. 



(4.4) 



(iii) 5*^(4, R) ^ S0{3, 2). The matrix A is a four-dimensional real matrix, satisfying 
A^ J = J A and TrA — 0, where 



J 



E 
-E 



E is the unit matrix of dimension 2. 

The bijection of an element (a;i, a;2, 0:3, 0:4, X5), of onto A is given by 



^ xi X2+X2, a;4 + 0:5 

X2 — —Xi —X4 — X5 

X4 — X5 Xl X2 — X3 

\ -X4 + X5 X2 + X3 -Xl 



(4.5) 



The transformation A' = SAS ^ with S G Sp{4, R) maps A into itself, namely, 
A'^J = J A', TrA' = 0. Since 

det A= {xl + xl-xl-xl+ xf) ^ = det^', (4.6) 

this transformation induces the desired isomorphism. 

(iv) Sp{l, 1) ^ 5'0(4, 1) A is defined by the four-dimensional complex matrix satisfying 
A'^J = JA, A^K = KA, Tr ^ = 0, with K = diag(l, -1,1,-1). 
The bijection of an element {xi,X2,X3,X4,X5) of R^ onto A is 



A 



Xl X2 + tXs 

X2 + iX3 —Xl 

X4 — ix^ 

\ —X4 + ix5 



X4 + ix5 \ 

—X4 — 1x5 

Xl —X2 + ixs 

X2 + 1X3 -Xl / 



(4.7) 



Given an element S of the group Sp{l, 1), that is to say, S^JS = J, S^KS = K, the 
transformation A' = SAS~^ maps A into itself. Since 



det A = — X2 — x^ - x^ — xlY = det A' 
this transformation induces the desired isomorphism. 



(4.8) 
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(v) SU{2,2) ^ SO{4:,2). A is defined by the four-dimensional complex matrix, 
satisfying = —A, A* I = I A, with A, the complex conjugate matrix of A, A* the 
dual matrix of A, namely, (A*)^^^ = ^SabcdA'^'^, and 



1 = 



E 
-E 



The bijection of an element {xi,X2, xs, X4, X5, xq) of onto A is (Beckers et al., 1978) 



A 



/ xi + 1x2 + 1x4^ X5 + ixQ 

— Xi — iX2 X^ — 1X6 —X3 + 1x4 

—Xs — iX4 —X5 + IXq —Xi + 1X2 

\ —X4 — ixQ Xs — 1x4 xi — 1X2 



(4.9) 



The transformation A' = SAS'^, with S satisfying S'^IS = /, maps A into itself. Since 
det A = [-xj -xl + xl + xl + xl + xf) = detA' (4.10) 

this transformation belongs also to 5*0(4,2). 

(vi) SL{2,Q) ^ 5*0(5, 1). Let A be a two-dimensional quaternion matrix defined by 
At = A. The biject ion of an element (xq, xi,X2,X2,X4,X5) of onto A is the following: 



A 



Xo + Xi X2+ Xsi + X4j + x^k 

X2 — Xsi — X4j — X^k Xq — Xi 



(4.11) 



with {i,j,k) a basis for the quaternions. The transformation A' = SAS"^, with S e 
SL{2, Q) maps A into itself. Since. 

det A= {xq-x\-xI-x\-x\- x\)^ = detA' (4-12) 

this transformation induces the desired isomorfhism (Barut et al., 1965) 

(vii) 5'L(4, R) ?s 5*0(3, 3). A is defined by the four-dimensional real matrix, satisfying 
A^ = —A. The bijection of an element {xi,X2,X3,X4, x^, xq) of R^ onto A is the following: 



A^ 



I —X\ -\-X4 X2 + X5 X3 + Xq \ 

—Xi — X4 X3 — Xq —X2 + X5 

—X2 — X5 —Xs + Xq Xi — X4 

\—Xs — Xq X2 — X5 —X1+X4 / 



(4.13) 



The transformation 5" = SAS'^, with 5* e 5'L(4, R), maps A into itself. Since 

det A= [xl+xl + xl- xl-xl- xf) = detA' 
this transformation induces the desired isomorphism. 



(4.14) 
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5. INTEGRAL TRANSFORMATIONS 



A transformation that leaves invariant an hypercubic lattice is called an integral 
transformation. It means that if we take an integral vector (the components of which 
are real integers or Gaussian numbers) the transformed vector is also integral. 

It is easy to prove that a transformation of a classical group is integral if and only if 
all its components are integers. 

Given the Cay ley realization of a classical groups (2.1) the corresponding transforma- 
tion is integral if all the Cayley parameters are integers and det {1 + X) = 1. But this 
procedure do not exhaust all the integral transformations. Let us give some examples. 



5.1 SO{S) 

The Cayley parametrization is given by (3.1.2.). Therefore we require m, n,p, q to be 
integer numbers and 

+ + = 1 (5.1.1) 

The non-exceptional matrix satisfying both conditions is the unit matrix correspond- 
ing to m = ±1, n = p = q = 

The other set of non trivial integral transformation are obtained in this way: from 
the Cayley transform 

it follows that S is integral if 2(1 + X)~^ is integral. 
From (3.1.2) we get: 

m? + q^ —mn + pq mp + nq 
mn + pq w? -\- p'^ —mq + np I (5.1.3) 



1 + X + + p2 + g2 



—mp + nq mq + np m^ + n^ 




Obviously 2(1 + X) ^ is integral if m, q are integers and m'^ + n'^ + p^ + q^ = 2. 
The solutions of this diophantine equation are: 

m — ±1 , n — ±1 , 

m = ±l , p=±l , n = q = 0} (5.1.4) 
m = ±1 , q = ±1 , 

By similar considerations we obtain also the integral transformations with 

m = ±l , n = ±l , p = ±l , q = ±l (5.1.5) 

5.2 ^0(3,1) 

As before we require in (3.3.2) all the Cayley parameters m, n,p, g, r, s, t to be integer 
numbers; after substituting A = {nt + ps + qr)/m the condition det (1 -|- X) — 1 becomes 

m^{m'^ + n^ +p^ + q^ - r"^ - -t^) - (nt + ps + qn)'^ = m^ (5.2.1) 
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This formidable diophantine equation can be simplified if we expand the inverse matrix 
appearing in the Caley transform (5.1.2), namely, 



= i-x + x^-x-^ + --- 

1 + X 

This series must be finite if matrix S is supposed to be an integral transformation, 
therefore = for some r. Since X is a nihilpotent matrix, it has zero as the only 
eigenvalue. Prom the secular equation it follows that the sum of all principal minors of the 
same order of the matrix X are equal to zero. This property applied to (3.3.1) gives the 
conditions 

(nt + ps + qr)'^ = y ■ ■ ) 

substituting these values in (5.2.1) we get m = 1. 

An other set of solutions are obtained from the condition 2(1 + X)~^ to be integral. 
This condition together with (5.2.2) gives 

2 2 

■X 



1 + X m? 

/ + — — —mn + pq + rs mp + nq + rt —mr = ns — pt \ 
mn + pq + rs + p^ — r"^ — t^ —mq + np + st —ms — nr + qt 

—mp + nq + rt mq + np + st m^ + n^ — r^ — s^ —mt + pr — qs 
\ —mr — ns + pt —ms + nr — qt —mt — pr + qs m^ + n'^ + p^ + q^ / 

(5.2.3) 

The condition of integral transformation requires that each matrix element multiplied 
by 2 must be divisible by m^. Therefore the only solutions are: 

m — 2 , n,t,p,s odd integers; q,r even integers 
or } (5.2.4) 

m = 2 , n,t,q,r odd integers; p, s even integers 

The same method can be applied to other integral transformations corresponding to 
classical groups. 



5.3. The method of isomorphism between real forms 

Section 4 can be used to find integral transformations. Take the isomorphism between 
SL{2,C) and S'0(3, 1). A spin transformation of the group SL{2,C) corresponds to a 
Lorentz tansformation with Cayley parametrization after the identification (3.3.4). 

An spin transformation is integral if the corresponding Lorentz transformation is 
integral. Schild has solved completely the problem of classifying all integral Lorentz 
tansformation by the following theorem (Schild 1949): 

A spin transformation (i,j = 1,2) is integral if and only if one of the following 
four conditions is satisfied: 

I. \ij are Gaussian integers such that A11A22 — A21A12 — 1 and such that An + A21 + 
A12 + A22 is even. 
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II. Xij = iJ-ij/il + i), where are odd integers such that yLtii//22 — A*2iA*i2 = 2i 

III. Xij are integers such that A11A22 — A21A12 = i and such that An + A12 + A21 + A22 
is even. 

IV. Xij = fJ'ij/il + i) where Hij are odd integers such that /X11//22 — /^i2A*2i = —2 
(For the definition and properties of Gaussian numbers see Schild 1949). 

We can compare this classification with our results of section 5.2. It is easy to 
check with the help of (3.3.4) that Case I corresponds to m = 1 and (5.2.2) and Case 
II corresponds to m = 2 and (5.2.4). 

Cases III an IV are obtained from I and II multiplying a, /3 or 7, 5 by i in (3.3.4). 
It corresponds to interchange first and second row in (3.3.2), multiplying the first one 
by -1. If we choose to multiply a, 7 or P,S by i in (3.3.4) we obtain the corresponding 
transformation (3.3.2) with the first and second column (multiplied by —1) interchanged. 

In order to complete the algorithm to calculate the SL{2, C) transformations with 
Gaussian numbers, we have to solve the diophantine equation A11A22 — A12A21 = 1. As 
before we write the general element of SL{2, C) as a finite series, say 

S = = 1-X with X^ = (5.3.1) 

1 + X ^ ^ 

the condition for X to be a nihilpotent matrix requires 

detX = trX ^0 , (5.3.2) 
the most general solution of which is 




(5.3.3) 



with z, w Gaussian integers. Therefore 

S^(l-zw, -zw'\^X^^ X^A 

\ Z, 1 + ZW J \X21 X22 J 

that can be enlarged to the general solution 

^^pn+pA2i Ai2+pA22\ (535) 
\ A21 A22 / 

with p, an arbitrary Gaussian integer. 



6. A LAGRANGIAN MODEL FOR THE KLEIN-GORDON FIELD ON THE 
LATTICE 

Let us define the real scalar field (pix, t) on the grid points of a (1 + 1) -dimensional 
lattice as (t>{sj, rn) = (f)^ where e, r are the fundamental space and time interval and j, n are 
integer numbers. We want to associate to this field a Lagrangian, such that the equations 
of motion are recovered from the Euler-Lagrange diff'erence equations. 
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A suitable Lagrangian for the Klein-Gordon field is: 



N-l 



Ln — 



N-l 

= eJ2^n (6.1) 

3=0 

where periodic boundary conditions for the field are supposed, Aj(Vj) are the forward 
(backward) difference operator and Aj(Vj) are the forward (backward) average operator 
with respect to the space indices: 



Aj/j = fj+i - fj ^jfj - fj - fj-i 



(6.2) 



A,/. = I {fj+i + h) V,/, = \ (/, + /,_i) 



(a;, t) 



'dt 



{x,t) 



(6.3) 



and similar for the time indices. In the limit j — > oo, e — > 0, js ^ x, n — > oo, r 
0, nr ^ t we have 



(6.4) 
(6.5) 

(6.6) 



Taking the variations of the Lagrangian density Cn with respect to the time difference 
of the field, we get 



dCr 



= -V„A,.V,>," 



VnA,.V,-7r; 



(6.7) 



with tt" as the conjugate momentum. 

To obtain the Euler-Lagrangian equation we take the time difference of the last 
expresion to be equal to the variation of the Lagrangian density with respect to the average 
field: 



1 



rAnVnA,Vj0^ = An 



1 



(6.8) 



where A„ has been introduce for homogeneity in the last equality and integration by parts 
have been used. 

The last expression is the wave equation for the Klein-Gordon field on the lattice. 

The "plane wave" solutions, the Fourier decomposition of the fields in terms of a 
complete set of orthogonal functions on the lattice have been given elsewhere (Lorente, 
1992). 
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With the help of the conjugate field defined in (6.7) we can construct an Hamiltonian 
density on the lattice in the usual way. 

Un = (VjAjVnTr^) ^VnAjVj4>] - (6.9) 
from which the Hamiltonian follows: 



if 2 1 



= eY^nn (6.10) 

i=o 

Taking the variation of the Hamiltonial density with respect to the averaged fields 
and its conjugate momentum, the Hamilton equations of motion are obtained: 

-An hiA^^n4>i] = V = (VjAjVnTT^) (6.11) 



dUn 



An \ ^VjAjVn4>] - M^VjAjVn4>] 



(6.12) 



where integration by parts has been used. Again, the Hamiltonian equations of motion 
lead to the wave equation (6.8). Notice that (6.11) and (6.12) can be simplified by the 
time average operator Vn, namely: 



^An (V,A^.0-) = An (v.A.-tt;) 
A. (V,-A,.7r;) = A. (^V,-A,.0- - M'VjA^cf>] 



(6.13) 
(6.14) 



which can be obtained from the following Hamiltonian 

Hn = eJ^l{ + il^.^-^i)' + M'{v,A,cl>]y\ (6.15) 

i=o ^ 

For the quantization of the Klein-Gordon field we introduce the equal time commu- 
tation relations: 



'3 J 



V,- A,-0-, VrAj,cj>} 1=0= V,- A,-7r7, V,. A,.^, 



(6.16) 
(6.17) 
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from which the Heisenberg equations of motion, the Fourier decomposition in terms of the 
creation and annihilation operators can be deduced in the usual way (Lorente 1992). 



7. A LAGRANGIAN FOR THE DIRAC FIELD ON THE LATTICE 

We can repeat the same steps as in the previous case for the Dirac field on the 
(1 + 1)— dimensional lattice {js, nr) = i^aj- 
A suitable Lagrangian density is 

Cn = -A,-A„V]" |747i ^A.-AnV," - ^^74A,•AnV'," + M^^AjAji;^"^ (7.1) 

with 

'0 -i\ fl \ . /O-l 



7i=(. , 74= _J, i7i74 = 75 = o' ^^-2) 



zA,-A^V^ = A.A^tt" (7.3) 



leading to the Euler-Lagrange equation 

d (iA^A.V',") 

-Vn i^AuT^tl = Vn .^^I" . = v4-A, A„Vj"747i " MA,A„V]"74| (7.4) 
^ ^ ^ a(A„A,V',"j U ' 'J 

Substituting (7.3) in (7.4), taking the adjoint operation of both sides and multiplying 
by 74 from the right we obtain the Dirac equation on the lattice 

Vn |7i^AjA„ - z74^AnAj + MA^Anj V" = (7.5) 

The plane wave solutions and the Fourier decomposition in terms of a complete set of 
orthogonal functions on the lattice have been given elsewhere (Lorente 1991). 

With the help of the conjugate field tt" = ■iV'j" we can construct the Hamiltonian 
density 

Hn = (AnA,-7r,") ^A„A,V; - 

= (a„A,V']") 7471 AnA.V'; + M (a„A,V']") Ar,A,^PJ (7.6) 

from which the Hamiltonian equations of motion can be derived leading again to the Dirac 
equation. 

As in the Klein-Gordon case, one can simplify the Hamilton equation by A„ which 
can be deduced from the new Hamiltonian 

N-l . . 

H^=eY, A,V']" -7471 A,V," + Mj.Aji;^ (7.7) 
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For the quantization of the Dirac field we require the equal time anticommutation 
relations 

[A,-C,-,A,.v£,]+ = J<^«pV (7.8) 

with other anticommutations vanishing. If we plagg in the Heisenberg equations of 
motion, we obtain the Dirac equation, from which the plane wave solutions and the Fourier 
decomposition in terms of the creation and annihilation operators can be obtained (Lorente 
1991). 



8. CONSERVATION LAWS AND LORENTZ IN VARIANCE 

As in the continuous case we can make the connection between symmetries and 
conservation laws in the language of generators. The condition for symmetry of the 
Lagrangian under space and time displacement and pure Lorentz transformation is that 
the generators are constant of the motion (Yamamoto 1991). In the case of the Klein- 
Gordon fields the generators of the (one step) space and time translations and Lorentz 
boost can be taken as: 

N-l 



(^^^-^7) iyAi^l) + (v,-4-^;) } (8.1) 



AT-l 

H = e 

j=o 

N-l 

K = e 

i=o 

N-l 



X; \ej I (a, V, A.tt;) ' + 1 (a, V,- A,0;) ' + (a,- V, A,.^^) ' 

7=0 

N-l -, 

- ^ E "^2 { (^^•^^■^.^) + (^.- ) } (8-3) 

Using (6.13) and (6.14) it can be proved that these operators are constant of time 

- A„P = -/^r^H = -AnK = (8.4) 

T T T 

In order to check the Lorentz invariance of quantized field scheme on the lattice, one 
can prove with the help of (6.16) and (6.17) that these operators satisfy the standard 
commutation relations: 

[H, P] = , [K, H]=iP , [K, P] = iH (8.5) 

For the Dirac quantum fields, the generators of the (one step) space and time 
translations and Lorentz boost can be taken as 

N-l 

^ = E (^^5) (^C,) (8.6) 
j=o 
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3=0 ^ 

N-1 N-1 



(8.7) 



(8.8) 



i=o i=o 
Using (7.5) and (7.8) it can be proved that these operators are constant of time 

-AnP = -AnH = -A^Mi4 = (8.9) 
r T T 

and that they satisfy the standard commutation relations 

[H,P]^0 , [H,Mi4]=P , [P,Mi4^H (8.10) 

We can convince ourselves that H and P are the generators of the time and space 
displacement by iteration of the Heisenberg equations of motion as in the continuous case. 
For the operator H we have 



A, (V,A,>;),if 



(8.11) 



where we have fixed the time index, say n = 0. 
By iteration we have 



t (:) (i)'(A„)' t (:) A?'(vA^?)„.„ 

fc=l ^ ^ ^ ' k=l ^ ^ 



k times 



which can be taken as the "Taylor expansion" on the lattice, namely. 



fc=i ^ ^ 



n=0 



(8.12) 



(8.13) 



In the limit n — > oo, r — > 0, nr ^ t the expression (8.12) becomes 

k=l ' ^ ^ ^ 7"^" ' fe=l 



k times 



= (j){x,t) 



t=0 

(8.14) 



But this expression is precisely the expansion of the continuous time traslations 
generated by the operator H 



e'"'(f){x, 0) e-'""' = (f){x,t) 
16 



(8.15) 



For the operator P we have 

1 r. 



where we have fixed the space index, say j = 0. By iteration we have 



(8.16) 



k-1 



fctimes 



k=l ^ ^ ^ 



which correspond to the j step space translation on the lattice. In the limit (8.17) becomes 
the continuous space translation generated by the operator P. 



APPENDIX. THE EINSTEIN DE BROGLIE RELATION ON THE LA- 
TTICE 

In order to make connection of our scheme with the Einstein-de Broglie relations 
E = hu>, p — hk we take the discrete plane waves solutions of (6.8). 

Obviously, we have for the period T and wave length A 

T = Nt , X = Ne {A.2) 

and for the phase velocity 

7 T 



If we impose the boundary conditions 

f^{k,u;) = f^{k,m) (A4) 
the wave number and the angular frequency can be defined as 

7 2 nm 2 nm ^ ^ at -i / ^ ,-\ 

km = -tg^, '^m = -tg^, m = 0,l,---,N -1 {A.5) 

Substituting the Einstein-de Broglie relations in the relativistic expresion E"^ — 
= (we use natural units ^ = c = 1), we obtain 
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Since the phase velocity and group velocity satisfy VpVg = 1, we have finally 
giving a discrete mass spectrum due to the lattice. 
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